Superradiance induced topological vortex phase in a Bose-Einstein condensate 
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We investigate theoretically a topological vortex phase transition induced by a superradiant phase 
transition in an atomic Bose-Einstein condensate driven by a Laguerre-Gaussian optical mode. We 
show that superradiant radiation can either carry zero angular momentum, or be in a rotating 
Laguerre-Gaussian mode with angular momentum. The conditions leading to these two regimes are 
determined in terms of the width for the pump laser and the condensate size for the limiting cases 
where the recoil energy is both much smaller and larger than the atomic interaction energy. 
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In addition to the spin angular momentum (SAM) as- 
sociated with its polarization [l[, electromagnetic radi- 
ation can also carry orbital angular momentum (OAM) 
associated with its spatial mode structure Q as for exam- 
ple in Laguerre-Gaussian (LG) beams It can exhibit 
quantum entanglement between OAM states and it 
was proposed [1] and demonstrated @ that superposi- 
tions of photonic OAM states can be utilized for higher- 
dimensional quantum communication and in dense data 
storage applications 0, Q . 

The efficient processing and storage of quantum infor- 
mation in terms of OAM requires reliable mechanisms 
for exchanging angular momentum between photons and 
atoms. For a supcrfluid state such as an atomic Bose- 
Einstein condensate (BEC), OAM states translate into 
topological excitations [91413, e.g. as vortices or circu- 
lating modes on a surface (l3|. The transfer of angu- 
lar momentum involving LG beam induced vortices in 
atomic condensates have been studied both experimen- 
tally [9T-fl2| and theoretically [H - 16 1 . The present paper 
reports conditions under which that transition can occur 
in a BEC simultaneously with the onset of a topological 
vortex. 

Under appropriate conditions, an ensemble of atoms 
optically driven above a threshold intensity can emit ra- 
diation in the form of a superradiant pulse (SR) [l?], [H| , 
a process analogous to a first-order phase transition 
(PT) [HEl. Early work on BEC superradiance 
analyzed it in terms of matter- wave gratings in the trans- 
lational [Hl-il [H, and in the polarization H|| 
degrees of freedom of the atoms. More recently sev- 
eral authors also discussed and demonstrated the onset 
of a structural phase transition from a homogenous to 
supersolid phase in a BEC trapped in an optical cav- 
ity [29|, (30[. In contrast, the present situation involves 
gratings in both translational and rotational degrees of 
freedom. Above the SR threshold, collective scattering of 
an incident LG pump laser results in the sudden and com- 
plete transfer of OAM to the BEC, bringing the conden- 



sate into a vortex state. Unlike the Raman coupled two- 
pulse pumping scheme of Refs. [l(J [H| this scheme yields 
several orders of magnitude larger vortex/no- vortex re- 
coil ratio. In addition, the recoiled atoms remain in the 



same internal state 11 



We consider a cigar-shaped BEC in a non-rotating 
elongated trap, pumped along the long condensate axis 
(z-axis) by a far off-resonant intense laser field of mo- 
mentum ko, see Fig. [TJ For a Fresnel number close to 
unity the atoms collectively recoil into well-defined mo- 
mentum states (side- modes), while the scattered light 
is predominantly along two end-fire modes propagating 
along ±z [II]. 

After adiabatic elimination of the excited states, the 
effective Hamiltonian describing the interaction of the 
condensate atoms of Bohr frequency cu a with a far-off 
resonant optical field detuned from the atomic transition 
by A = uj a - uj is 

H = [ d 3 r^(r)H g (r)^(r)+J2 I <2 3 k/ia; fc aj^a M 



d 3 rd 3 kd 3 k'g u ,(k, k'; r)^(r)4,a k ^(r)(,l) 



where H g (r) is the atomic Hamiltonian, ip(r) and the an- 
nihilation operator for atoms in their electronic ground 
state, and fik,£ are optical field mode annihilation oper- 
ators, the indices £,£' = 0, ±1 labeling the angular mo- 
mentum of the optical modes 0] . 
The effective coupling coefficients 



<(k,k';r) = - 



hg*(k)g(k<) 
A 



§J,,(r)*k',<'(r) (2) 



are determined by the single atom-photon dipole matrix 
element g(k). Here, $k,f (r) are the mode functions of the 
light field with wave number k and angular momentum 
M. For LG modes we have 
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of width we which carries £K units of OAM along z-axis. 

It is known that supcrradiance in an atomic BEC can 
be modeled by an effective Hamiltonian that only in- 
cludes the dominant end-fire modes of the light field and 
the associated matter- wave side- modes [15|, |26|, |33j . We 
can then proceed by expanding the condensate field op- 
erator in terms of side modes 



26| as 



(4) 



where y m , 9z (r) = Lp m {r)e" n ^e lqzZ is the eigenfunction 
of an atom with recoil momentum q z in the z-direction 
and vortex charge m. It is given by the solution of the 
equation H m {r)ip m {r) = fl m ip m (r) where 

fr , ^ ( d2 1 d \ m2fl2 rrt N 

H ^ V) =2M {^ + rTr) + 2M^ + Vt{r) > (5) 

and Vt is the trap potential. For a harmonic trap, these 
states are similar to the LG modes 11511 
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where w is the radial width of the condensate. 

Focusing on first-order scattering processes reduces 
then the Hamiltonian of the system to the simplified form 



m,q z I J 

/ dk '9im(kz,qz)cl tim al rA a^\co t o + h.c.{7) 



where e m (q z ) = h(tt m + w R (q z )), is the LG pump 

laser mode, Co,o describes the initial condensate, and 

/, \ t -9*(kz)g(k ) 

g£7n(k z ,q z ) = —n — <V+m,l 

x J d 3 r^ )m (r)$^^(r)$ feo)1 (r)^o,o(r). (8) 

When driven by a Gaussian beam of wave vector ko, 
supcrradiant scattering occurs mainly in two counter- 
propagating end-fire modes of wave vectors k ~ zfckn, 
where ko is the wave vector of the pump photon. (22l.[23|. 
with the initial condensate coupled to two dominant 
matter- wave side modes of momenta q = and q = 2ko. 
With a LG driving field, in contrast, four end- fire opti- 
cal modes are excited, with k ~ ±ko and I = or 1, 
accompanied by four matter-wave side modes k ~ ±ko 
with m = 1 or 0, the optical end-fire modes with £ = 
being coupled to matter-wave modes with m = 1 due to 
conservation of angular momentum (see Fig. [1]). 

The dynamics of scattering from a supcrfluid BEC fun- 
damentally differs in the two regimes. If the recoil en- 
ergy (jvjjn{q) = h 2 q 2 /2M) is much smaller than the in- 
teraction energy (per atom), single particle excitations 



(a) 

C = -i<o . . 

•4st/v/%, <T q = 2k„ Ji\t 

t = l m=0 


(b) 

(=0 m =l 


(c) 

iT,=+k„ 


(d) 

k,=+k„ 


m=0 ' = 1 





FIG. 1: (Color online) Cigar-shaped BEC initially in the state 
m = 0, q = and illuminated by a strong LG-mode laser of 
wave vector ftko and in the OAM mode 1=1 and. Superradi- 
ant scattering is predominantly along the end-fire modes with 
wave vector k e = ±ko. Four possible end- fire modes and the 
corresponding matter- wave side-modes are shown in (a)-(d). 



are highly suppressed [llj, |31|, l32| . The collectivity, in 



duced by the interparticle interaction, does not allow 
atoms to recoil one by one, i.e. induction of vortex only 
after a critical rotation frequency 13[ in rotating trap. 
In this mean-field regime, photons are scattered by only 
the quasiparticle excitations and BEC is described by an 
order parameter ^(r). For relatively large recoil energy 
[32| , condensed atoms behave as non-interacting ones and 
photon scattering is essentially due to individual atoms 



one by one [lfj, [111, |2lH23 



In the large recoil energy regime, rotatory forward scat- 
tering (Fig. Id) of LG laser is neglected on the basis 
that (i) it cannot transfer sufficient linear momentum 
and excitation energy }34L l35| and (ii) such low-energy 
scatterin g ar e suppressed by the structure factor of the 
BEC 3l|,|36|. In the mean-field (smaller recoil frequency) 



regime, where single atom scattering is already forbidden 
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361 ] . all processes Fig. l(a-d) included. 
Since the large recoil energy regime corresponds to a 
series of well-known experiments j2l!423| , we investigate 
it first. Keeping only relevant modes results in the effec- 
tive Hamiltonian (h = 1) 

^int = - (5l4fco,O a -fco,i a fco i ,l £o - 0+h - C -) 

- (s24fco,i a -fco 1 O a Sl«D.O + h - C -) +? ^> ( 9 ) 

where c 2 k ,o, c 2 fe 0) i, a-fc ,i and a_ fcoj0 are the matter- 
wave side modes and optical end-fire modes illustrated in 
Figs.Q]} andQJi. Here gi = gio = w 2 w 2 /2it 2 (w 2 + w 2 ) 2 
and 92 = .901 = {wi/w)g\ can be calculated from Eq. ([8]). 

and 7id = — 2gocJ o4 ' "H^O * s a diagonal light shift 
term with go = gi- The free-field terms have been elim- 
inated through a second rotating frame transformation 
after moving to a co-rotating frame at the laser fre- 
quency [13]. The first term in Eq. (|9]) yields normal 
(non-rotatory) SR, while the second term induces vor- 
tex excitations in the BEC and results in rotatory SR. 
For times short enough that the depletion of the ini- 
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tial condensate can be ignored we can make the substi- 
tution co.o v^V, where N is the number of condensed 
atoms. In that limit the side-mode population dynam- 
ics can be treated analytically [2^], giving dc2k ,m/dt ~ 
G m Nc2k , m /2, with G m oc g^+v F «r g 2 /gi = W\jw > 1, 
the initial exponential growth of vortex side-mode occu- 
pation from initial fluctuations is faster than that of the 
non-vortex side mode. This suppression of normal SR 
relative to rotatory SR, in the large recoil energy regime, 
can be thought of as a topological vortex phase transi- 
tion. For example, a typical ratio 52/91 = 2 fioj ] would 
result in vortex/no- vortex side mode population ratio of 
3 x 10 4 for TV = 10 6 . We note that the condition w\ > w , 
is consistent with OAM transfer into classical objects [37[ 
and with the two-pump Raman vortex excitation proto- 
col in BEC [lfj. In these cases, however, the resulting 
vortex/no- vortex population ratio increases linearly with 
(wi/w) 2 , in contrast to the exponential growth charac- 
teristic of the present situation. 

To investigate the small recoil energy regime, we in- 
clude the effects of two-body collisions and we invoke the 
mean-field approximation whereby the field operators are 
replaced by c- numbers, specifically ip — >• ip, a^k ,i ot\, 

S-feo,o «2, afco,o ~> "3, and a^i ~* «l, without side- 
mode expansion ([?]). With these approximations, the 
Heisenberg equations of motion reduce to 
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Here Uq = <7i/A, g$ = ^Kha s /m, a s being the s-wave 
scattering length, — Ai^^ are the end-fire mode frequen- 
cies in the rotating frame at frequency uq, and 



a, 13 



d 3 r^ akojm {r)^ ko<m ,(v)\^(r,t)\ 2 , (14) 
±1 label the sign of the wave vectors 



where a and (3 - 
of amplitude fco . 

Figures [2] and [3] summarize the results of a typical 
numerical solution of these equations. In this example 
Wi = 1.2w, NU Q = 10, r\ = a L VNU = 15, Ai = 1, 
A2 = 1.001, and A3 = 10~ 3 in units of the recoil fre- 
quency wr. These values are comparable to those of 
Ref. |2l|, |30(. (We note that the energy scale of light 
shift is several orders of magnitude larger than all other 
terms, however its effect is transient in and of no signifi- 
cant importance for the time scales under consideration. 

) 
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FIG. 2: (Color online) Snapshots illustrating the onset of BEC 
SR in a BEC driven by a LG laser beam, (a) The BEC 
starts in a Gaussian state; (b) Onset of normal SR, with the 
condensate density mapping the LG intensity profile, (c) Just 
before the onset of rotatory SR the axial (^-symmetry of the 
condensate is broken, with a higher density near the <j> ~ n 
angle, (d) Rotatory SR induces a vortex topological phase 
transition in the BEC. 



The results of our simulations are perhaps most dra- 
matically illustrated in Fig. [21 which shows four snap- 
shots of the transition from an initial condensate in the 
ground state of the harmonic trap to a condensate with a 
density profile that maps the intensity profile of the LG 
driving field, but still with zero average angular momen- 
tum, and subsequently to a rotatory condensate carrying 
one quantum of angular momentum per atom, that is, a 
condensate in a vortex state. 

The dynamics of the transition arc shown in more de- 
tail in Fig. [3J which plots the evolution of key observ- 
ables during the SR transition to the vortex phase: For 
the parameters of these simulations the normal SR PT 
happens at t ~ 8/ojr (Fig- [2h)- It coincides with a peak 
in the intensity |ai| 2 (Fig. [3^ of the LG end-fire mode 
fl-Aio.i- Spatial order sets in along the z-direction, sec 
(Fig. [3}j), as evidenced by the non-zero value of (e 2 ** 02 ) 
and the fact that the condensate acquires linear momen- 
tum {p z )/hkQ = 2|ai| 2 as a result of momentum conser- 
vation. In this early phase, there is no OAM transfer to 
the condensate happens as both the laser and the end-fire 
modes have equal winding number £ = I. 

This early dynamics is followed by the onset of rota- 
tory SR. As a precursor to that transition the conden- 
sate first exhibits an axial symmetry breaking both in its 
phase and density, as illustrated in Fig. This is the 
first indication of OAM transfer from the optical field to 
the atoms [HI . Rotatory SR reveals itself with a sudden 
increase in the intensity |o!3 1 2 of the Gaussian end-fire 
mode a_fc .o, shown in Fig. [3^, and a decrease in the 
Laguerre-Gaussian mode |o;i| 2 . At this point the angu- 
lar momentum is transferred fully from the optical field 
to the matter wave, establishing a vortex in the BEC 
(Fig. EH), the BEC gaining (L z ) = Nh (Fig. [3fc) of or- 
bital angular momentum. Comparison of Figs.[3lt and [3b 
indicates that the total OAM of the system is conserved 
in this transition, during which the process the linear mo- 
mentum along z is transferred back to the optical field, 
as can be seen from Figs. [3^ and[3j3. 

We finally remark that we found numerically that the 
critical pump rates (ry) are different for the onset of non- 
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FIG. 3: Time evolution of the condensate for the parameters 
of Fig. [2] (a) Mean photon number in the Gaussian end-fire 
mode |ct3 1 2 and in the LG end-fire mode, (b) Linear momen- 
tum transfer ((p 2 ) oc |ai| 2 ) and expectation value of the two- 
photon momentum recoil operator {e 2lk ° z }; (c), Expectation 
value of the normalized angular momentum {L z )/Nfi. 



rotatory ( m = 9.5) [H H , and rotatory SR {i] 2 = [20 
12.4): Below rji no SR scattering occurs. Between the 
two values (771 < ?/ < 772) only normal SR takes place, 
and rotatory SR occurs for above 772- 

Summarizing, we have studied theoretically the mu- 
tual induction of a topological vortex phase transition 
in a BEC and a rotatory superradiant phase transition- 
like . The cooperative nature of SR from an incident 
LG laser beam allows for a sudden transfer of a large 
amount of OAM into a condensate, bringing it into a 
vortex state. For most current condensate SR experi- 
ments, non-rotatory SR can be exponentially suppressed 
by rotatory SR if the transverse width of the LG pump 
laser is larger than the condensate transverse width. 
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